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Abstract 

We study the ^5*0 superfluidity of A hyperons in neutron star matter and neutron stars. We use 
the relativistic mean field (RMF) theory to calculate the properties of neutron star matter. In the 
RMF approach, the meson-hyperon couplings are constrained by reasonable hyperon potentials that 
include the updated information from recent developments in hypernuclear physics. To examine the 
^5*0 pairing gap of A hyperons, we employ several AA interactions based on the Nijmegen models 
and used in double- A hypernuclei studies. It is found that the maximal pairing gap obtained is a 
few tenths of a MeV. The magnitude and the density region of the pairing gap are dependent on 
the AA interaction and the treatment of neutron star matter. We calculate neutron star properties 
and find that whether the ^5*0 superfluidity of A hyperons exists in the core of neutron stars mainly 
depends on the A A interaction used. 

PACS numbers: 26.60.-c, 21.30.Fe, 24.10.Jv 

Keywords: Superfluidity, Relativistic mean field theory. Neutron star matter 



Electronic address: songtc@nankai . edu. on 



I. INTRODUCTION 



Neutron stars are natural laboratories for studying the physics of dense matter. They 
are associated with some of the most exotic environments in the universe. Our knowledge of 
neutron star interiors is still uncertain. The central density of neutron stars can be extremely 
high, and many possibilities for such dense matter have been suggested fl-3|. For densities 
below twice normal nuclear matter density (po ~ 0.15 fm^'^), the matter consists of only 
nucleons and leptons. When the density is higher than 2po, the equation of state (EOS) 
and composition of matter are much less certain. The presence of hyperons in neutron stars 
has been studied by many authors |4-^. K~ condensation in dense matter was suggested 
by Kaplan and Nelson [9] and has been extensively discussed in many works ^ y,llO|. It 
has been suggested that the quark matter may exist in the core of massive neutron stars, 
and the hadron-quark phase transition can proceed through a mixed phase of hadronic and 



,ua.k .a«e. |ll a H If deconfined oua.k ,na«e. does ex.=t in.de =ta., . likely to 
be in a color superconducting phase pj, 16[ , and various color superconducting phases have 



been intensively investigated in recent years 



12 



-ll4|. 



Baryon pairing is believed to play an important role in the evolution of neutron stars 
HQ. The presence of neutron superfluidity in the crust and the inner part of neutron stars 
can be considered well established [15!]. The neutron fluid in the crust probably forms a ^5*0 
superfluid. With increasing density, the ^Sq interaction turns repulsive, and the neutrons 
in the outer core mainly form a superfluid. On the other hand, one can expect ^5*0 
proton pairing in the outer core because the small proton fraction brings about a low proton 
density in this region. In the inner core of neutron stars, hyperons may appear through the 
weak interaction because of the fast rise of the baryon chemical potentials with density. It 
is widely accepted that hyperons appear around 2po- The presence of hyperons tends to 
soften the EOS at high density and lower the maximum mass of neutron stars 4M8ll as well 

nn 

as increase the neutron star cooling rate p^, 117[. In general, the flrst hyperon to appear is 
A, which is the lightest one with an attractive potential in nuclear matter. The potential 
of S hyperons is now considered to be repulsive; therefore S~ appears at a higher density 
than A in neutron star matter [l8, Q. The ^Sq superfluidity of A hyperons is suggested to 
occur in the same way as that of neutrons arising from the attractive AA interaction in the 
^5*0 channel 20l-l24|. It is known that hyperon pairing can signiflcantly affect the thermal 
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evolution of neutron stars by suppressing neutrino emission from the hyperon direct Urea 



process 



24j-|26l| . Youn g ne utron stars cool primarily by neutrino emission from the interior. 



As discussed in Refs. 



27 



29|, the neutrino emissivity in superfluid matter is exponentially 



suppressed when the temperature T is much lower than the superfluid critical temperature 
Tc. On the other hand, superfluidity initiates a specific neutrino emission from the Cooper 
pair breaking and formation process, which is forbidden in nonsuperfluid matter. This 
process is exponentially suppressed when T <^ Tc, and it is much less efficient than the 



direct Urea process [29|, ISOj. Hence the presence of baryon superfluidity can drastically 



suppress the neutrino emission, which may play a key role in neutron star cooling. We 
are mainly interested in the possibility of ^5*0 superfluidity of A hyperons in neutron stars. 
So far, the ^5*0 pairing gap of A hyperons is still uncertain because it can be significantly 
influenced by both the properties of matter and the AA interaction. More studies are needed 
to determine these uncertain factors using available information from recent developments 
in hypernuclear physics. 

In this article, we focus on the ^5*0 superfluidity of A hyperons in neutron star matter, 
which is composed of a chemically equilibrated and charge-neutral mixture of nucleons, 
hyperons, and leptons. To calculate the pairing gap, we need to specify how to treat the 
neutron star matter and the AA interaction. In this article, we use the relativistic mean field 
(RMF) theory to calculate the properties of neutron star matter. The RMF theor y h as been 
successfully and widely used for the description of nuclear matter and finite nuclei 31n35|. It 
has also been applied to providing the EOS of dense matter for use in supernovae and neutron 
stars {sgJ. In the RMF approach, baryons interact through the exchange of scalar and 
vector mesons. The meson-nucleon coupling constants are generally determined by fitting 
to some nuclear matter properties or ground-state properties of finite nuclei. To examine the 



influence of the RMF parameters, we employ two successful parameter s ets, TMl [37[ and 



NL3 [38^, which have been widely used in many studies of nuclear physics 



As for 



the meson-hyperon couplings, there are large uncertainties because of limited experimental 
data in hypernuclear physics. Generally, one can use the coupling constants derived from the 
quark model or the values constrained by reasonable hyperon potentials. The meson-hyperon 
couplings play an important role in determining the properties of neutron star matter |8l. |40|. 
We use the values constrained by reasonable hyperon potentials that include the updated 
information from recent developments in hypernuclear physics. We take into account the 
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two additional hidden-strangeness mesons, a* and 0, which were originally introduced to 



obtain the strong attractive AA interaction deduced from the earlier measurement 41|. 



42 



A recent observation of the double-A hypernucleus AAHe, called the Nagara event 
has had a significant impact on strangeness nuclear physics. The Nagara event provides 
unambiguous identification of ^A^e production with a precise AA binding energy value 
-Baa = 7.25 ± 0.19]'^o;}^ MeV, which suggests that the effective AA interaction should be 
considerably weaker (Ai?AA — 1 MeV) than that deduced from the earlier measurement 
(ASaa — 4-5 MeV). The weak hyperon-hyperon (YY) interaction suggested by the Nagara 
event has been used to reinvestigate the properties of multistrange systems, and it has been 
found that the change of YY interactions affects the properties of strange hadronic matter 



dramatically 



45| . We would like to examine whether the ^5*0 superfluidity of A 



hyperons exists in neutron star matter, and how large the pairing gap can be if it does, by 
considering recent developments in hypernuclear physics. 

The aim of this article is to investigate the possibility of forming a A superfiuid in neu- 
tron stars. It has been suggested that hyperon superfluidity could significantly suppress the 
neutrino emission in the core of a neutron star and play a key role in neutron star cool- 



ing 



24 



25| . Over the last decade, there has been some discussion in the literature about 



hyperon pairing in dense matter 



20 



24 



46| . In the work of Balberg and Barnea 20|], the ^5*0 



superfluidity of A hyperons has been studied by using an effective AA interaction based on 
a G matrix calculation and an approximation of nonrelativistic effective mass obtained from 
single-particle energies. Their calculation predicts a gap e nergy of a few tenths of a MeV for 
A Fermi momenta up to about 1.3 fm~^. In Refs. 2l|, |22|, Takatsuka and Tamagaki studied 
A superfluidity using two types of bare AA interactions based on the one-boson-exchange 
(OBE) model and two types of hyperon core models. They found that A superfluidity could 
exist in a density region between 2po and (2.6-4.6)po, depending on the pairing interaction 
and the hyperon core model. A study of AA pairing in a pure neutron background has been 
presented by Tanigawa et al. 23| using the relativistic Hartree-Bogoliubov model, where the 
AA pairing gap was found to decrease with increasing background density and decreasing 
AA attraction. Both A and superfluidities in neutron star matter have been investigated 
by Takatsuka et al. {24 1 using three pairing interactions based on the OBE model and several 
nonrelativistic EOS with different incompressibilities. It was found that both A and S~ are 
superfiuid as soon as they begin to appear at around 4po, although the pairing gap and the 
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density region depend on the pairing interaction and the EOS of neutron star matter. The 



effect of the Nagara event on A superfluidity has been discussed in Refs. [23|, |2J] , where the 
weak attractive AA interaction suggested by the Nagara event leads to very small pairing 
gap in dense neutron matter {23] or the disappearance of A superfluidity in neutron star 



matter 



24| . All these studies indicate that the AA pairing gap in dense matter depends 



both on the AA interaction and on the EOS of matter. The key role of A superfluidity in 
neutron star cooling motivates us to investigate the possibility of forming a A superfluid in 
neutron stars by carefully considering the pairing interaction and the description of neutron 
star matter with the updated information from recent developments in hypernuclear physics. 

This article is arranged as follows. In Sec. [ITl we briefly describe the RMF theory for 
the calculation of neutron star matter properties. In Sec. IIIIl we discuss the AA interaction 
used in the gap equation. We present the numerical results in Sec. |IVl Section |V] is devoted 
to a summary. 



II. RELATIVISTIC MEAN FIELD THEORY 



We use the RMF theory to describe the neutron star matter, which is composed of a 
chemically equilibrated and charge-neutral mixture of nucleons, hyperons, and leptons. In 
the RMF approach, baryons interact through the exchange of scalar and vector mesons. The 
baryons considered in this work are nucleons {p and n) and hyperons (A, S, and S). The 
exchanged mesons include isoscalar scalar and vector mesons {a and u), an isovector vector 
meson (p), and two additional hidden-strangeness mesons {a* and (p). The total Lagrangian 
density of neutron star matter takes the form 

Crmf = [hfid'' -rriB- QaBO - ga*BCr* - QcoBlfiUj'' 

B 

-g^Blt,4>^ - QpSlt^TiBpt] i^B + ^d^ad^'a - ^mla^ 
+ ^C3 {u^uj^' f - ^Rif.yR'^'' + ^mjpi^pf 

^^i}i{%-i^,d^ -mi\i)u (1) 
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where V'b and ipi are the baryon and lepton fields, respectively. The index B runs over the 
baryon octet {p, n, A, S"*", S^, and the sum on / is over electrons and muons 

(e~ and fi^). The field tensors of the vector mesons, u, p, and 0, are denoted by W^^^, Ri^u, 
and 5*^,^, respectively. In the RMF approach, the meson fields are treated as classical fields, 
and the field operators are replaced by their expectation values. The meson field equations 
in uniform matter have the following form: 

ml. + 9.a- + 9sa' = - E ^ ^ .Tl ..^"^^^ (2) 
m^u + C3U =2_^ — , (3) 



B 

37r2 



B 

^2 * = _y9^ —^^k'dk, (5) 

where a = (a), oj = (oj^), p = (pso), o"* = (o"*), and = (0°) are the nonvanishing 
expectation values of meson fields in uniform matter; m*^ = uib + (^o-bct + ga*BO'* is the 
effective mass of the baryon species B, and kp is the corresponding Fermi momentum. 

The meson-baryon coupling constants play an important role in determining the proper- 
ties of neutron star matter. To examine the infiuence of the RMF parameters, we employ 



two successful parameter sets, TMl 



371 and NL3 



in the present calculation. These 



parameters have been determined by fitting to some ground-state properties of finite nuclei, 
and they can provide a good description of nuclear matter and finite nuclei, including un- 
stable nuclei. With the TMl (NL3) parameter set, the nuclear matter saturation density is 
0.145 fm~^ (0.148 fm~^), the energy per nucleon is —16.3 MeV (—16.3 MeV), the symm etry 
energy is 36.9 MeV (37.4 MeV), and the incompressibility is 281 MeV (272 MeV) 



As for the meson-hyperon couplings, we take the naive quark model values for the vector 
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coupling constants; 

1 1 1 

gpN = -gpT, = gpE, gpA = o, 

2fi'</,A = 25f</,E = g<f,E = — —guiN, g<pN = 0. (7) 

The scalar coupling constants are chosen to give reasonable hyperon potentials. We de- 
note the potential depth of the hyperon species i in the matter of the baryon species j by 
Uji''\ It is estimated from the experimental data of single-A hypernuclei that the potential 



depth of a A in saturated nuclear matter should be around ^ —30 MeV 



47j. For 



S hyperons, the analysis of S atomic experimental data suggests that S-nucleus potentials 
have a repulsion inside the nuclear surface and an attraction outside the nucleus with a 
sizable absorption. In recent theoretical works, the S potential in saturated nuclear mat- 
ter is considered to be repulsive with a strength of about 30 MeV l^^]. Some recent 
developments in hypernuclear physics suggest that H hyperons in saturated nuclear matter 



have an attraction of around 15 MeV 



18|, 



48J. In this article, we adopt f/j^^ = -30 MeV, 



Uj^ = +30 MeV, and f/4 = —15 MeV to determine the scalar coupling constants. We 
obtain, for the TMl (NL3) parameter set, g^A = 6.228 (6.323), g^^ = 4.472 (4.709), and 
fi'o-s = 3.114 (3.161), respectively. The hyperon couplings to the hidden-strangeness meson 



a* are restricted by the relation U^^ ^ uP ^ 2Ui^^ ^ 2Uj{'^ obtained in Ref. |49fl. The 



weak YY interaction implied by the Nagara event suggests uj^^ ~ —5 MeV, and hence we 
obtain g^,^ = 5.499 (5.678) and g^*s = 11-655 (11.899) for the TMl (NL3) parameter set. 
We assume ga*j] = ga*A and take rria* = 980 MeV and = 1020 MeV in this article. 

For neutron star matter consisting of a neutral mixture of baryons and leptons, the /3- 
equilibrium conditions without trapped neutrinos are given by 

l^p = Ais+ = /in - /ie, (8) 

/^A = /^E" = /iH'J = /in, (9) 
/is- = /iH- = /in + /ie, (10) 
fJ'P = /ie, (11) 

where /ij is the chemical potential of species i. At zero temperature, the chemical potentials 
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of baryons and leptons are given by 

fJ'B = \J kp'^ + m*j^ + gujBOJ + g^B4> + gpBTsBP, (12) 
Hi = k^p^ + , (13) 

respectively. The electric charge neutrality condition is expressed by 

Pp + Ps+ = Pe + + Pe- + Ph- , (14) 

where pi = (kp)^ /(3vr^) is the number density of species i. We solve the coupled Eqs. ([2])- 
©, fl8l)- f|TIl) . and f|T^ self-consistently at a given baryon density ps- Then we can calculate 
the EOS and the composition of neutron star matter as well as the effective mass and the 
Fermi momentum of A hyperons, which are crucial in the study of AA pairing. 



III. GAP EQUATION AND AA INTERACTION 

We study the ^Sq superfluidity of A hyperons in neutron star matter. The crucial quantity 
in determining the onset of superfluidity is the energy gap function A{k), which can be 
obtained by solving the gap equation 

. (15) 

^[E{k')-E{k^)f + A^k') 

where E {k) is the single-particle energy of A with momentum k. For A hyperons in neutron 
star matter, the single-particle energy in the RMF approach is given by 



E{k) = sj k'' + mf + g^KOJ + 

The effective mass m\ and the Fermi momentum kp are computed self-consistently at a 
given baryon density ps within the RMF approach. 

For the A A pairing interaction in the ^5*0 channel, the potential matrix element can be 
written as 

V {k,k') = {k\V^^{^S,)\k') = A-K J r'dr Mkr)VAA{r)jo{k'r), (16) 

where jo{kr) = sin{kr) / (kr) is the spherical Bessel function of order zero and Vaa (^) is 
the ^5*0 AA interaction potential in coordinate space. Because of large uncertainties in the 
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AA interaction, we adopt several AA potentials. Most of them are based on the Nijmegen 
models and are used in double-A hypernuclei studies, which are of the three-range Gaussian 
form 

3 

^AA(r) = 5^i;.exp(-rV/3,2). (17) 

i=l 

The short-range term provides for a strong soft-core repulsion, whereas the medium-range 
and long-range terms provide for attraction. The parameters Vi and /3j are taken from 



Refs. 



50 



53| . and we list them in Table [H The NDl potential was given in Ref. 50| as 



an effective soft-core interaction fitted to the Nijmegen model D (ND) hard-core interac- 
tion. Another simulation of the ND interaction, called ND2 in this article, was obtained 
in Ref. |5l|. The ESCOO, NSC97b, NSC97e, and NSC97f potentials given in Ref. [sj were 
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53| were obtained by adjusting 



obtained by changing the strength of the medium-range attractive component of the three- 
range Gaussian potential such that they could reproduce the scattering length and the 
effective range as close to values by the corresponding Nijmegen models. It is well known 
that the Nagara event provides unambiguous identification of aa^g production with a pre- 
cise AA binding energy value i?AA and has had a significant impact on strangeness nuclear 
physics. The NFs and NSC97s potentials given in Refs. 
parameters to reproduce the experimental value of i?AA(AAHe) from the Nagara event. We 
have also chosen an Urbana-type potential that has been successfully used to ex pla in the 
experimental values of hypernuclei. The Urbana potential could be found in Ref. 54 1. 

We plot in Fig. [T] all AA potentials considered in the present work. The strongest AA 
interaction is the ESCOO potential, whereas the weakest AA interaction is the NSC97f po- 
tential. We note that the NFs, NSC97s, and Urbana potentials simulate the experimental 
value of -BaaIaaHg) from the Nagara event. 



IV. RESULTS AND DISCUSSION 

In this section, we investigate the ^5*0 superfluidity of A hyperons in neutron star matter 
and neutron stars. We employ the RMF model with the parameter sets TMl and NL3 
to calculate the properties of neutron star matter, which is known to provide excellent 
descriptions of the ground states of finite nuclei, including unstable nuclei. The meson- 
hyperon couplings play an important role in determining the properties of neutron star 
matter. We use the values constrained by reasonable hyperon potentials that include the 
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updated information from recent developments in hypernuclear physics. As for the AA 
pairing interaction used in the gap equation, we adopt several AA potentials that have been 
used in double- A hypernuclei studies. Some simulate the experimental value of -BAA(AAHe) 
from the Nagara event. With the effective mass and the Fermi momentum of A hyperons 
obtained in the RMF approach, the gap equation [Eq. (IT^ J is solved numerically. 

In Fig. [21 we show the resulting ^5*0 pairing gap of A hyperons at the Fermi surface, Ap, 
as a function of the baryon density, pB, in neutron star matter. The results of TMl and 
NL3 are plotted in Fig. [2] (top) and Fig. [2] (bottom), respectively. In the case of TMl (NL3), 
the threshold density of A is around 0.31 fm~^ (0.28 fm~^), and A hyperons form a ^5*0 
superfiuid as soon as they appear in neutron star matter. With increasing baryon density, 
increases first, reaching a maximum value at pb ~ 0.34 fm~^ (p^ ~ 0.30 fm~^), then 
decreases and finally vanishes at pB < 0.46 fm~'^ (p^ < 0.38 fm~^) for the case of TMl 
(NL3) with the ESCOO potential. It is found that the maximal pairing gap is about 0.8 
MeV with the ESCOO potential in the TMl case. This is because the ESCOO potential has 
the strongest attraction among the AA interactions used here. The pairing gaps with the 
NDl and ND2 potentials are of the order of 0.1-0.2 MeV, as shown in Fig. [21 In addition, 
we find that the pairing gaps are of the order of 10"^ MeV (TMl) or absent (NL3) with 
the NSC97e, NFs, NSC97s, and Urbana potentials. The A pairing does not appear for the 
NSC97b and NSC97f potentials. We present in Table [TTl the maximal pairing gap at the 
Fermi surface (A™*") and the corresponding baryon density (pb), effective A mass {m\), 
and Fermi momentum (kp) using these potentials with the TMl and NL3 parameter sets. 

The A pairing gap Ap depends not only on the AA interaction but also on the properties 
of A hyperons in neutron star matter. In Fig. [3l we show the particle fraction, Yi = Pi/pB, as 
a function of the baryon density, pB, using the RMF model with the TMl (Fig. [3l top) and 
NL3 (Fig. [31 bottom) parameter sets. It is seen that A hyperons appear around 0.31 fm~^ 
(TMl) or 0.28 fm~^ (NL3) and then increase rapidly with increasing density. We note 
that hyperon threshold densities, fractions, and effective masses are dependent on the RMF 
parameters used. This dependence has an effect on the resulting pairing gap, as shown 
in Fig. [2l Our results with the NDl potential can be compared with those in Table III 



of Ref . 



24| . where the same AA interaction (called the ND-Soft potential) was used. The 



difference is in the treatment of neutron star matter, for which they use a nonrelativistic G 
matrix-based effective interaction approach, whereas we use the RMF approach. In our case 
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of TMl (NL3), the maximal pairing gap at the Fermi surface is 0.17 MeV (0.12 MeV), as 
given in Table [TTl where ps = 0.344 fm'^ (ps = 0.303 fm"^), Ya = 0.039 (Fa = 0.044), and 
m\ = 743 MeV {m\ = 706 MeV). Takatsuka et al. 2J] obtained the maximal pairing gap 
of 0.34 MeV at ps = 4.5po ioT the TNI6u EOS. The larger pairing gap at higher or given in 



Ref. 



2J] is because of their smaller Ya and larger m^. As discussed in Refs. [20|, l23|, |2J] , the 



pairing gap is very sensitive to the effective mass. Generally, a smaller effective mass leads 
to a higher single-particle energy and then yields a smaller pairing gap. In Fig. |U we show 
the effective mass of A hyperons, function of the baryon density, ps, using the 

RMF model with the TMl (solid line) and NL3 (dashed line) parameter sets. It is shown 
that decreases with increasing ps- When A hyperons appear around 0.31 fm^^ (TMl) 
or 0.28 fm~^ (NL3), the effective mass of A hyperons is about 762 MeV or 727 MeV. It is 
found that the effective masses of A hyperons in the NL3 case are smaller than those of 
TMl, which leads to a smaller pairing gap, as shown in Fig. [2l We note that the effective 
mass of A hyperons is mainly determined by the coupling constants g^A and ga*A- Here we 
use the values constrained by reasonable hyperon potentials, which are suggested by the 
experimental data of single- A hypernuclei and by the Nagara event. 

To examine whether the ^5*0 superfluidity of A hyperons exists in neutron stars, we solve 
the Tolman-Oppenheimer-Volkoff (TOV) equation with the EOS of the RMF theory over a 
wide density range. For the nonuniform matter at low density, which exists in the inner and 
outer crusts of neutron stars, we adopt a relativistic EOS based on the RMF theory with 
a local density approximation 0, The nonuniform matter is modeled to be composed 
of a lattice of spherical nuclei immersed in an electron gas with or without free neutrons 
dripping out of nuclei. The low-density EOS is matched to the EOS of uniform matter 
at the density where they have equal pressure. The neutron star properties are mainly 
determined by the EOS at high density. Using the EOS described in Sec. [TTl we calculate 
the neutron star properties and find that the maximum mass of neutron stars is about 1.70 
Mq (2.06 Mq) with the TMl (NL3) parameter set. According to the compilation of measured 
neutron star masses js, 55|, some massive neutron stars were reported to be observed recently. 
However, the uncertainties in these mass measurements are rather large, and the mass of PSR 
J0751 + 1807 was corrected from (2.1 ± 0.2) Mq to (1.26 ± 0.14) M© [5Q]. We note that the 
EOS used here could not be ruled out by current observations. In Figs. |5] and |6l we show the 
central baryon density as a function of the neutron star mass. We find that whether the ^5*0 
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superfluidity of A hyperons exists in the core of neutron stars depends on the AA interaction 
used. With weaker AA interactions, such as NSC97b and NSC97f, the A superfluidity does 
not appear inside neutron stars. For the NSC97e, NFs, NSC97s, and Urbana interactions, 
although we obtain the pairing gaps of the order of 10^^ MeV in the TMl case, it is unhkely 
that A superfluidity can exist in observed neutron stars because of its low superfluid critical 
temperature ~ 0.57Air/KB ~ 10^ K 



24 



29|. With stronger AA interactions, such as 



ESCOO, NDl, and ND2, the ^5*0 superfluidity of A hyperons may exist in massive neutron 
stars, as shown in Figs. |5] and [61 In the case of TMl (NL3) with the ESCOO potential, A 
hyperons do not appear in neutron stars with M < 1.37 (M < 1.50 Mq). For neutron 
stars with 1.37 Mq < M < 1.63 Mq (1.50 Mq < M < 1.82 Mq), A hyperons in the core of 
neutron stars form a ""^^o superfluid. However, when M > 1.63 Mq (M > 1.82 Mq), not only 
superfluid A but also normal (nonsuperfluid) A can exist in the core of neutron stars because 
the central baryon density exceeds the upper limit of the range where A superfluidity exists. 
The presence of nonsuperfluid A hyperons in the core of massive stars would lead to a more 
rapid cooling than the case with only superfluid A hyperons. The mass region, where only 
superfluid A hyperons exist in the core of neutron stars, is shaded in Figs. |5] and El It 
is shown that the region with the ESCOO potential is the widest among all cases in these 
flgures. This is because the ESCOO potential has the strongest attraction, and its pairing 
gap covers the widest density range, as shown in Fig. |2l We note that this region depends 
both on the AA interaction and on the EOS of neutron star matter. 



V. SUMMARY 

We have studied the ^Sq superfluidity of A hyperons in neutron star matter and neutron 
stars. In this article, we employ the RMF model with the parameter sets TMl and NL3 
to calculate the properties of neutron star matter, which is composed of a chemically equi- 
librated and charge-neutral mixture of nucleons, hyperons, and leptons. The RMF theory 
has been successfully and widely used for the description of nuclear matter and flnite nuclei, 
including unstable nuclei. In the RMF approach, baryons interact through the exchange 
of scalar and vector mesons. The baryons considered in this article are nucleons {p and 
n) and hyperons (A, E, and H). The exchanged mesons include isoscalar scalar and vector 
mesons (cr and u), an isovector vector meson (p), and two additional hidden-strangeness 
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mesons {a* and 0). It is well known that the meson-hyperon couplings play an important 
role in determining the properties of neutron star matter. We have used the couplings con- 
strained by reasonable hyperon potentials that include the updated information from recent 
developments in hypernuclear physics. To examine the ^5*0 pairing of A hyperons, we have 
adopted several AA potentials. Most are based on the Nijmegcn models and have been used 
in double-A hypernuclei studies. NFs, NSC97s, and Urbana potentials have simulated the 
experimental value of B\A{\jJie) from the Nagara event. 

We have calculated the ^-S'o pairing gap of A hyperons at the Fermi surface, Ap, using 
the AA potentials adopted in this article. It is found that depends both on the AA 
interaction and on the treatment of neutron star matter. The maximal Ap obtained in the 
present calculation is about 0.8 MeV with the ESCOO potential in the TMl case. This is 
because the ESCOO potential has the strongest attraction among the AA interactions used in 
this article. The NDl and ND2 potentials yield somewhat smaller Ap of the order of 0.1-0.2 
MeV. For the NSC97e, NFs, NSC97s, and Urbana potentials, the values of are of the 
order of 10"^ MeV (TMl) or absent (NL3). The A pairing does not appear for the NSC97b 
and NSC97f potentials. The difference in these results reflects the dependence of Aj? on the 
AA interaction. On the other hand, the magnitude and the threshold density oi Ap are also 
dependent on properties of neutron star matter, especially on the effective mass and particle 
fraction of A hyperons. In the case of TMl (NL3) with the ESCOO potential, the threshold 
density of Ap is around 0.31 fm~^ (0.28 fm~"^), reaches a maximum value at ps ~ 0.34 fm~^ 
{Pb ~ 0.30 fm~^), and finally vanishes at ps < 0.46 fm~^ {pB < 0.38 fm~^). By solving 
the TOV equation, we have calculated neutron star properties and found that whether the 
^5*0 superfluidity of A hyperons exists in the core of neutron stars mainly depends on the 
AA interaction used. With stronger AA interactions, such as ESCOO, NDl, and ND2, the A 
superfluidity may exist in massive neutron stars. It is unlikely that A superfluidity can exist 
in neutron stars with the NFs, NSC97s, and Urbana interactions, which have simulated the 
experimental value of B\A{%jJie) from the Nagara event. 

In this article, we have considered the updated information from recent developments in 
hypernuclear physics and used the weak attractive AA interactions suggested by the Nagara 
event. However, there are still large uncertainties in the hyperon- hyperon interaction and 
the EOS of neutron star matter. A more precise study of the A pairing in neutron stars 
requires further development in hypernuclear physics. 
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TABLE I: Parameters of ^Sq AA interaction defined in Eq. (jl7p . taken from Refs. 5[H53|. The 
size parameters are the same for all cases, which are /3i = 1.342 fm, /32 = 0.777 fm, and /^a = 0.350 
fm. The strength parameters are in MeV. 





Vi 


V2 


^"3 


NDl 


-21.92 


-283.5 


4745 


ND2 


-21.49 


-379.1 


9324 


ESCOO 


-21.49 


-456.6 


9324 


NSC97b 


-21.49 


-182.1 


9324 


NSC97e 


-21.49 


-207.1 


9324 


NSC97f 


-21.49 


-177.1 


9324 


NFs 


-10.96 


-141.8 


2137 


NSC97S 


-21.49 


-250.1 


9324 
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TABLE II: Maximal pairing gap at the Fermi surface A™^^ obtained with several AA potentials; 
Pb the total baryon density of neutron star matter where A™^^ is obtained, and kp and are 
the corresponding Fermi momentum and effective mass of A hyperons, respectively. 







TMl 










NL3 






PB 


i,A 




A max 




PB 


fvp 


ml A^^=^ 




(fm-3) 


(fm-1) 


(MeV) 


(MeV) 




(fm-3) 


(fm-1) 


(MeV) (MeV) 


NDl 


0.344 


0.738 


743 


0.17 




0.303 


0.731 


706 0.12 


ND2 


0.339 


0.681 


747 


0.10 




0.298 


0.664 


711 0.06 


ESCOO 


0.349 


0.789 


740 


0.81 




0.305 


0.762 


704 0.62 


NSC97b 


















NSC97e 


0.329 


0.548 


753 


1.2 X 10" 


-4 








NSC97f 


















NFs 


0.329 


0.548 


753 


5.4 X 10" 


-4 








NSC97S 


0.329 


0.548 


753 


4.0 X 10- 


-4 








Urbana 


0.329 


0.548 


753 


5.5 X 10" 


-4 










0.6 0.8 1.0 1.2 1.4 1.6 
r (fm) 

FIG. 1: (Color online) ^Sq AA interaction potentials used in this article. 
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FIG. 2: (Color online) pairing gap of A hyperons at the Fermi surface A^t^ as a function of 
baryon density ps in neutron star matter with the NDl, ND2, and ESCOO potentials: (top) TMl 
and (bottom) NL3. 




FIG. 3: (Color online) Particle fraction Yi = Pi/ pb as a function of baryon density pB'- (top) TMl 
and (bottom) NL3. 
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FIG. 4: (Color online) Effective mass of A hyperons m\ as a function of baryon density ps- 
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FIG. 5: Central baryon density pc as a function of neutron star mass M in the TMl case. The 
region where only superfiuid A hyperons exist in the core of neutron stars is shaded. 
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FIG. 6: Same as Fig. [U but for the NL3 case. 
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